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Abstract. A non-perturbative local moment approach to single-particle dynamics 
of the general asymmetric Anderson impurity model is developed. The approach 
encompasses all energy scales and interaction strengths. It captures thereby strong 
coupling Kondo behaviour, including the resultant universal scaling behaviour of the 
single-particle spectrum; as well as the mixed valent and essentially perturbative 
empty orbital regimes. The underlying approach is physically transparent and innately 
simple, and as such is capable of practical extension to lattice-based models within the 
framework of dynamical mean-field theory. 



1. Introduction 

The Anderson impurity model [1] (AIM) has long played a pivotal role in understanding 
the physical behaviour of materials dominated by strong, local Coulomb interactions (for 
a comprehensive review see e.g. [2]). A broad range of problems is encompassed by the 
AIM itself, including [2] magnetic impurities in metals, heavy fermion systems when 
coherence effects are suppressed, and the burgeoning area of quantum dots [3,4]. 

Further impetus to the study of quantum impurity models has arisen with the 
advent of dynamical mean-field theory [5-8] (DMFT), within which correlated lattice- 
based systems such as the periodic Anderson or Hubbard models, reduce to an effective 
quantum impurity hybridizing self-consistently with the surrounding fermionic bath. To 
capture such problems entails the ability to describe an AIM with essentially arbitrary 
dynamics (w-dependence) in the hybridization function A(a>) that embodies coupling 
between the impurity and the underlying host/bath. The theoretical difficulties here 
are considerable, unsurprisingly given the wide range of physics encompassed; and 
despite impressive progress in recent years [5-8] there remains a need for new theoretical 
approaches that can handle in particular the full range of interaction strengths: from 
non-perturbative strongly correlated regimes with their attendant low-energy scales, 
through to weak coupling, essentially perturbative domains of behaviour. 

Similar comments are in fact applicable even to the conventional metallic AIM 
[1]. Here, static properties (thermodynamic and related) are certainly well understood 
using a variety of powerful approaches, including the numerical renormalization group 
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[9], Fermi liquid theory [10] and the Bethe ansatz [11]. But the situation is less 
satisfactory when it comes to a theoretical description of dynamics, in particular single- 
particle excitations. A wide variety of theories, approximate of necessity, have of 
course been developed in this regard [2]; including for example perturbation theory in 
the interaction strength (U) [12], self-consistent renormalization thereof [13], modified 
perturbation theory [14], large- N expansions [15,16], the non-crossing approximation 
[17-19] and generalizations of it to finite-?/ [20,21], slave boson mean- field theory 
[22,23], a conserving t-matrix approximation [24] and the spinon approximation [25,26]. 
Their undoubted virtues notwithstanding however, many of these approaches have well 
recognized qualitative limitations [2]; and there is much scope for further theoretical 
development, ideally via an approach that is sufficiently general and practicable that it 
can be extended with relative ease to handle correlated lattice-based models within the 
DMFT framework. 

We have recently initiated development of one such, the local moment approach 
(LMA) [27-32]. This non-perturbative method is technically simple and transparent, 
with the physically intuitive notion of local moments [1] introduced explicitly and self- 
consistently from the outset. For the symmetric AIMs in which context it has thus far 
been considered, the LMA handles single-particle dynamics for all interaction strengths, 
and on all energy scales including recovery of Fermi liquid behaviour at low energies. It 
captures in particular the spin-fluctuation physics symptomatic of the strong coupling 
Kondo regime, manifest in an exponentially narrow Kondo resonance in the single- 
particle spectrum D(uj) [27]; such that the resultant scaling behaviour D(u) = F(uj/ook) 
(with uk the Kondo scale) can be obtained in closed form [28], and gives excellent 
agreement with NRG calculations [33] that provide benchmark results against which 
to compare approximate theories. The LMA has also been extended to finite-T [31], 
encompassing both single-particle dynamics and associated transport properties such as 
the resistivity. The role of a magnetic field H [30], which poses particular difficulties 
for conventional theoretical approaches to dynamics, can likewise be addressed; and 
for the Kondo regime in particular, corresponding static properties [29] such as the 
Wilson ratio, impurity magnetization and spin susceptibility, are found to agree well 
for essentially all field strengths with exact results known e.g. from the Bethe ansatz. 
Finally but importantly, we add that the LMA is not restricted to the Fermi liquid 
physics arising ubiquitously in the metallic AIM. The soft-gap AIM [34], containing an 
underlying quantum phase transition between generalized Fermi liquid and degenerate 
local moment phases, provides a particular example; in which a rich range of behaviour 
has been uncovered by the LMA [32] and confirmed by NRG calculations [33]. 

The LMA considered hitherto is nonetheless specific to particle-hole symmetric 
AIMs where the impurity orbital energy, e;, is slaved to the on-site interaction via 
€i — It is obviously desirable to extend the approach to encompass the generic 

asymmetric case, where the enlarged parameter space spans a wider range of physical 
behaviour [2], from the strong coupling Kondo domain through mixed valent behaviour 
to the ultimately perturbative empty orbital regime. This is important also from the 
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viewpoint of extension to lattice-based models within DMFT [5-8] , where asymmetry in 
the underlying impurity model corresponds to lattice-models away from half-filling; and 
is hence required to capture e.g. heavy fermion physics employing the periodic Anderson 
model. 

A local moment approach to the asymmetric AIM is considered in the present 
paper. Following the requisite background (§2), the two-self-energy description inherent 
to the LMA is considered in §3, with particular emphasis on the notion of (self- 
consistent) symmetry restoration that is central to the approach. That discussion is 
general, being applicable to an essentially arbitrary diagrammatic approximation for 
the underlying self-energies, and not dependent upon the symmetry-specific arguments 
hitherto employed for the symmetric AIM [27,28]. The particular non-perturbative 
approximation to the LMA self-energies implemented here in practice is discussed in 
§4.1; it passes the criterion of practicability, yet appears to capture rather well the 
relevant regimes of behaviour. Results arising therefrom are given in §5, including the 
rather subtle issue of universal spectral scaling in the strong coupling Kondo limit (§5.1); 
as well as evolution of single-particle spectra from strong to weak coupling behaviour 
(§5.2) as the Kondo, mixed valence and empty orbital regimes are traversed. The paper 
concludes with a brief summary/outlook. 



2. Background 

With the Fermi level taken as the energy origin, the familiar Hamiltonian for the spin-| 
AIM [1,2] is 

H = ^2 e ><nk* + J^(ei + \Uhi- a )h ic + V lk {c\ a c ku + cJ, CT c i(T ). (2.1) 

k,a a k,a 

The first term describes the non-interacting host with dispersion e k , and the third is the 
one-electron host-impurity coupling. The second term refers to the correlated impurity 
with site-energy e x and on-site interaction U; and the large parameter space of the model 
is conveniently specified by the asymmetry parameter [1] 

^=l + f (2-2) 

with i] = for the particle-hole symmetric case (for which e\ = — ^, and the impurity 
charge n = J2 a (n ia ) = 1 for all U). 

Our interest is in single-particle dynamics embodied in the impurity Green function 
G{ijj){^ G(t) = —i(T(c ia (t)c\ a ))), which is naturally independent of spin a (=| / j 
or +/— ) since H is invariant under a <-> —a; and hence the impurity spectrum 
D(uj) = — ^sgn(cj)ImG'(c<j). In the trivial non-interacting limit U = 0, the Green 
function reduces to 

g(u)=[u+-e i -A(u)}- 1 (2.3) 
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where u + = u + iO + sgn(co>), and is determined solely by the hybridization A (a;) 
A R (u) - isgn(w)A I (a;) given by 



A M = E^r^- ( 2 - 4 ) 

The essence of the conventional AIM is that the host is metallic by presumption, 
corresponding to a non-zero hybridization strength defined by A = Aj(u = 0) (with 
uj = the Fermi level). In practice, and without any essential limitation, we thus 
consider the usual wide flat-band host [2] for which A(u;) = — isgn(u;)Ao; where the 
hybridization strength is related by A = nV 2 p to the host density of states p and 
the matrix elements V = VJfc. The model is thus characterized by two independent 
parameters: either 

*-£ 6 - x ( 2 - 5 > 

^0 ^0 

as is traditional; or, as may prove more convenient (and in fact essential in the Kondo 
scaling regime as we show later), by fixed asymmetry 77 and one or other of U , t\. 
For U > 0, G(uj) is conventionally written as 

G{u)=[g{u)- 1 -V{u)Y x (2.6) 

in terms of the single self-energy S(cu) = S r (c<j) — isgn(a;)S I (u;), with E R /E I related by 
Hilbert transformation. The limiting \ow-u behaviour of G(u) that is symptomatic of 
the Fermi liquid character of the AIM, is the familiar quasiparticle form [2] 

= UJ+/Z -e[ + isgnH [A + 0(u 2 )} (2 ' 7) 

that is simply the leading low-u expansion of equation ( ^.6[ ) using ~ 0(uj 2 ); 

where Z = [l — (<9E R (u>) / du) w= o\ is the quasiparticle weight, and e- = e\ + S R (0) the 
renormalized level. The latter is related to the excess charge n imp induced by addition 
of the impurity, via the Friedel sum rule [35] (which itself follows directly from the 
Luttinger integral theorem [36]); specifically [2] 

n imp = 2 f du Ap(u) = 1 - -tan" 1 ( ^±IlM\ ( 2 .8) 

such that (e- =) e\ + S R (0) = A tan[|(l — n imp )]. Here Ap(u) = p{uj) — phost(^) is the 
change in total density of states of the system due to addition of the impurity, given by 
Ap(uj) = — ^sgn(a;)Im {G(lj) [1 — dA(u)/duj]}; so that n imp is related in general to the 
local impurity charge 

n = 2 I du D(u) (2.9a) 
J —00 

by 

n imp = n - -Im [ du G(u)^^- (2.9b) 
7T du 
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(with n imp = n for the wide flat-band case considered in practice). Equations ( |2.7| ,8) 
lead directly in turn to 

nA D(uj = 0) = sin 2 (fn imp ) (2.10) 

relating the Fermi level spectrum to n imp , and reducing to the 'pinning' condition 
itAqD(0) = 1 for the particle-hole symmetric case where n; mp = 1 always. The above 
results, all well known, will prove important in subsequent sections. 

The quasiparticle form equation (|2.7|) , conjoined with the Friedel sum rule for 
e[ = €i + S R (0), forms a starting point for microscopic Fermi liquid theory as well 
as being the essential result for G(u) arising from slave boson mean-field theory (see 
e.g. [2]). Its familiarity should not however obscure its limitations, for it is confined 
to the lowest frequencies u/A Z <C 1; and as such (save trivially for U = 0) captures 
only a tiny fraction [28] of the Abrikosov-Suhl resonance, let alone high-energy spectral 
features such as the Hubbard satellites. To capture dynamics on all energy scales, 
while recovering correctly the limiting quasiparticle form equation ( |2.7|) , is not a trivial 
matter; and non-perturbative approaches are certainly required to handle e.g. the strong 
coupling Kondo regime. One such is discussed here. 

Before proceeding, two remaining relevant points should be made. First, briefly, 
to specify the minimal range of asymmetry i] (equation ( |2.2j )) that need be considered. 
Under a particle-hole transformation, e y is replaced by — [e; + U] [9]; and, labelling 
temporarily the ej-dependence of G{uj) and n imp , it is straightforward to show that 
G(u;e i ) = -G(-u;-[e i + U)),i.e. 

D(u>;e{) = D(-u;-[e i + U]) (2.11) 

and hence n imp (e;) = ni mp (— [e; + U]). From this it follows that only e; > —U/2 (and 
hence < n imp < 1) need be considered, corresponding (equation Q2.2Q ) to r\ > 0. This 
range is assumed from now on. 

The second point refers to the strong coupling Kondo limit of the AIM. 
Corresponding to n imp — > 1, this arises when e ; < such that |ei| ^> 1 and U — |ei| ^> 1 
(i.e. the lower/upper Hubbard satellites centred on e x = — and e; + £7 respectively are 
well below/above the Fermi level). The approach to the Kondo limit is not therefore 
unique, in that n imp — > 1 arises for any given asymmetry 77 = 1 — 2\ii\/U G [0, 1] upon 
progressively increasing either U or |e ; |. This is reflected in turn by the fact that the 
Kondo model onto which the low-energy sector of the AIM maps in strong coupling 
under a Schrieffer- Wolff transformation [37], contains both exchange (J) and potential 
(K) scattering contributions [2]; viz 

H K = e kfi ka + 2 Js ; • S(0) + l -K c L c fc'- ( 2 - 12 ) 

fc,CT fc,fc',cr 

Here S; and S(0) denote respectively the impurity spin and conduction electron spin 
density at the impurity; and the exchange and potential scattering matrix elements are 
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related to the bare parameters of the AIM by [2]: 



Ao < 




7T 


M\ 


Ao < 




7T 





u-h\ 



(2.13a) 



pK = ^{ — - w -j- l j. (2.13b) 

The appropriate Kondo model may thus be characterized by the two dimensionless 
parameters pJ and Kj J, with p J <C 1 for the AIM — > Kondo reduction to be valid; and 
where from equations (|2.13a| ,b) and ( ^2[ ) 



t = 1 -it^ < 2 - 14 > 

is simply the asymmetry of the underlying AIM. Note that the Kondo model that in 
practice is usually considered lacks potential scattering, and hence corresponds to the 
strong coupling limit of the symmetric AIM alone (77 = 0). In general however, physical 
properties of the AIM in the strong coupling regime should depend on Kj J — rj; 
a statement that, with one limiting exception, applies in particular to the scaling 
behaviour of the single-particle spectrum, D{uj) = F(u/uk) with u K oc A Z — > the 
Kondo scale. That exception resides in the leading low-u; behaviour of G(u) embodied 
in the quasiparticle form G(ou), equation ( |2.7|) ; which in the Kondo limit where n imp — > 1 
(and hence e- = Cj + S R (0) — > 0), is dependent solely upon u/A Z oc uj/uk with no 
explicit dependence on Kj J = rj. To our knowledge one nonetheless has no a priori 
reason to expect that the single-particle scaling spectrum on all energy scales uj/ujk will 
be independent of the asymmetry 77. 

3. LMA: basis 

The conventional route to single-particle dynamics is via the usual 'single' self-energy 



£(u;). But this is merely defined by the Dyson equation implicit in equation (|2.6| ), and 
a determination of G(u) in this way is not obligatory; indeed there may be good reasons 
to avoid it, notably the inability of conventional perturbation theory to handle strong 
correlations in general. The LMA thus eschews such an approach, focusing instead on 
a two-self-energy description that is a natural consequence of the mean-field approach 
from which it starts. In this section we consider briefly the implications of such, in a 
manner that is neither dependent on the essential details of practical implementation 
(§4.1) nor confined to the particle-hole symmetric case considered hitherto [27-32]. 

There are three essential elements to the LMA. (i) First that local moments ('/i'), 
viewed as the primary effect of interactions, are introduced explicitly from the outset. As 
in Anderson's original work [1], the starting point is thus simple static mean-field (MF), 
i.e. unrestricted Hartree Fock. This contains in general two degenerate broken symmetry 
MF states (reflecting the invariance of H under a <-> — a); denoted by a = A or B and 
corresponding respectively to local moments p = +\p\ and — \p\. Notwithstanding 
the severe limitations of MF by itself, it may nonetheless be used as a basis for a 
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genuine many-body approach encompassing the correlated electron dynamics that are 
the essence e.g. of Kondo physics, (ii) To this end the LMA employs the two-self-energy 
description that follows naturally from the underlying two mean-field saddle points; with 
non-trivial dynamics introduced in practice (§4.1) into the associated self-energies via 
coupling of single-particle excitations to low-energy transverse spin fluctuations, (iii) 
The final, central notion behind the LMA is that of symmetry restoration: self-consistent 
restoration of the broken symmetry endemic at mean-field level, and recovery of Fermi 
liquid behaviour, as pursued below. 

Within the LMA, G(u) is expressed formally as 



a 



(3.1a) 



with propagators G a(J {ijj) = [g -1 (u;) — S^u;)] -1 built from the appropriate MF 
state a = A or B, and self-energies separated as Tj arj {uj) = E° CT + £ QCT (u;) into a 
purely static contribution TP aa that alone is retained at pure MF level; together with 
Tj a(T {ijj) = Sco-^Qo-}] that in particular contains the key dynamics, and is a functional 
of (and built diagrammatically from) the underlying MF propagators Qaa{^>)- The 
first, brief issue here is rotational invariance: the fact that G(u) is independent of 
spin, er. This is correctly preserved, since the invariance of H under o —a implies 
Ga<t(u) = Gb-o-(^); whence the sum in equation ( |3.1a| ) is indeed a- independent. By 
the same token, equation (|3.1a|) may be written as 



G{u) 



(3.1b) 



involving a 'spin sum' that is independent of a. Equations ( 3.1a| ,b) are entirely 
equivalent. We choose to work with the latter form and, since the a label is redundant, 
we drop it from now on. 

The impurity Green function is thus 

G{u)= l -[G ] {u) + G i {u)} (3.2a) 



with 



G„ 



-1 



and self-energies T> a (uj) = T^(u) — isgn(a;)E^.(a;) separated as 



(3.2b) 



(3.3) 



with S° given by the static bubble diagram (and 'everything else' in E CT (u;)). 

The conventional single self-energy E(a;) follows directly from an underlying two- 
self-energy description, being given on comparison of equations ( |2.6| ) and (|3.2a| ,b) by 

- i [i(£ t (w)-Ei(w))r 
E T (w) + + 



9 



-i 



;o;)-|(E T (a;) + S i (a;)) 



(3.4) 
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with g{ui) the non-interacting propagator equation ( |2.3| ). At the pure MF level of 
unrestricted Hartree Fock (HF), the dynamics of Ti a (uj) are neglected and £ CT (u;) = S° = 
y(n — cr|/i|) (with the local charge (n) and moment determined self-consistently in 
the usual way, §3.1); from equation ( |3.4| ) the corresponding self-energy is then: 

v ( \ U (^l/^l) 2 

Shf(^) = -n + — r , y i r^- 3.5 

2 cj+ - |_e; + ^-nj - A(u;) 

The problem with MF by itself is clear, for if \fi\ ^ then H l (ou = 0) ^ and Fermi 
liquid behaviour is violated. This is not surprising, for the resultant degenerate local 
moment state is not perturbatively connected to the non-interacting limit; but it is no 
less wrong for that and, despite the physical appeal of local moment formation, is a 
major reason why MF has not hitherto provided a successful basis for a many-body 
approach. This problem does not of course occur if = is enforced (restricted 
HF), but there another one arises. The two- and single- self energy descriptions then 
coincide, and S H f(^) = ^ n (= ^°) is just the static Hartree contribution. This produces 
a trivial energy shift to the non-interacting propagator, and subsequent construction of 
the dynamical £(w) via conventional perturbation theory in U is essentially equivalent 
to expanding about the restricted HF state. But when local moments can form at 
MF level, this saddle point — in contrast to those of unrestricted HF — is unstable 
(a local maximum, with the single restricted HF determinant unstable to particle- 
hole excitations). And it is this that underlies at heart the divergences that plague 
conventional perturbation theory [2] if one attempts to perform the more or less standard 
diagrammatic resummations (e.g. RPA) that one might expect are required to capture 
the regime of strong correlations. We do not of course doubt the applicability in principle 
of perturbation theory in U — for the metallic AIM; but practice is another matter. 

The LMA in a sense seeks the best of both worlds: to retain the two-self-energy 
description with the notion of local moments (and essential stability of the underlying 
MF state), while incorporating dynamics into the associated self-energies T, a (uj) in such 
a way that the resultant description is simple and tractable, and yet recovers Fermi 
liquid behaviour at low energies. 

3.1. Symmetry restoration 

The first question here is the latter point: under what conditions on the 'E a (u) , s will 
the single £(u;) exhibit Fermi liquid behaviour at low-cu, i.e. E I (o;) ~ 0{u 2 )l To answer 
this, consider a simple low-u; expansion of the E CT (a;)'s (equation ( |3.3|) ), viz 



'■l-i 
z„ 



uj (3.6a) 



for the real parts, where Z a — [1 — (<9E^(co>) /doj)^^]" 1 is thus defined and no a priori 
constraints are imposed on (0) = E° + T, a (uj = 0) at the Fermi level uj = 0; together 
with 

= ) ~ a^ 2 (3.6b) 
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for the imaginary part. The latter form is guaranteed from the diagrams (see §4) 
for S -(c<j) = So-fj^o-}] provided the host is metallic (A 7^ 0, as appropriate to the 
metallic AIM), but we emphasize is not sufficient by itself to guarantee ~ 0(u 2 ). 

Equations ( |3.6a| ,b), together (if appropriate) with a low-o; expansion of the hybridization 
A(o>), may then be used in equation (|3.4j ) to determine the Iow-cj behaviour of the 
resultant single self-energy. This is a matter of algebra, and we find thereby that the 
necessary/sufficient condition for 

~ 0(oj 2 ) (3.7) 

is that: 

tf(oo = 0) = E R (w = 0). (3.8) 

If equation (|3.8| ) is satisfied, then from equations ( p.4|) and ( |3.6a| ,b) (i) all self-energies 
coincide at the Fermi level, 

E r (uj = 0) = E R (cu = 0) (3.9) 

(for either a); (ii) the low-u; behaviour of E r (c<j) is 

E R (cu) ~ S R (0) - - lj u (3.10) 

with the usual quasiparticle weight Z = [l — (dY, R (u) / 800)^=0] 1 related to the {Z a } 
by Z^ 1 = \{Z7 + Z^ 1 ); and (iii) the quasiparticle form equation ( |2.7[ ) for G(u) is 
recovered. 

Equation ( |3.8j ), referring solely to the Fermi level u = 0, is the symmetry restoration 
(SR) condition that is central to the LMA; and which, if satisfied, guarantees Fermi 
liquid behaviour. It is quite general, meaning not confined to the particle-hole symmetric 
AIM considered hitherto (which is recovered as a particular case of the above, but 
was originally argued for [27,28] on different, symmetry-specific grounds). The general 
consequences of SR in practice are nonetheless as found for the symmetric model [27,28] : 
Self-consistent imposition of equation ( |3.8|) amounts (§4) to a self-consistency condition 
for the local moment \n\ and, most importantly, generates a low-energy spin-flip scale 
u m . The latter, manifest in particular as a strong resonance in the transverse spin 
polarization propagator lmll + ~(uj), is the Kondo scale, exponentially small in strong 
coupling; its physical significance in the approach being that it sets the timescale 
t ~ h/u m for restoration of the broken symmetry inherent at crude MF level (and 
arising in effect from dynamical tunneling between the degenerate MF minima). We 
also add that if the SR condition equation ( |3.8| ) cannot be satisfied, then a doubly 
degenerate local moment phase results [32] (with a characteristic spin-flip scale u m = 
reflecting the local degeneracy, and hence r = 00). While this does not arise for the 
metallic AIM where SR is ubiquitously satisfied, it is the self-consistent possibility of 
such embodied in SR that enables the LMA to capture [32] e.g. the quantum phase 
transition from a (generalized) Fermi liquid to a local moment state in the soft-gap 
AIM (Ai(u) ~ |cL>| r ) [34] where both phases arise. 
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Finally, note that with the SR condition equation 
rule equation (|2.8|) may be written as 



e ; + = 0) = A tan 



-(1 



n 



imp J 



satisfied, the Friedel sum 



(3.11) 



(independent of a). And for later use, the LMA G CT (u;)'s may then be expressed using 
equations ( |3.2b| ) and ( |3.3| ) as 



G„ 



OJ 



— A(u) — A tan 



-(1 



n 



imp ) 



[E a {u) - E CT (0)) 



-1 



(3.12a) 



being dependent only on the 'non-MF' contributions (E^u;)) to H a (u); such that n imp 
is given via equation (2.9) (using equation (3.2)) by 



n 



imp 



7T 



GJu) 



1 - 



dA(u) 
du 



(3.12b) 



The Friedel sum rule is not however satisfied by approximate theories in general, 
save for the particle-hole symmetric case where it is guaranteed by symmetry. That 
is true e.g. of the NCA [38], and even for conventional second order perturbation 
theory in U [12], where n imp inferred from e; + E R (0) = A tan[|(l — n imp )] (equation 
(p.8|)) using the second-order E R (0) does not coincide with n imp obtained from spectral 
integration. Satisfaction of the Friedel sum rule, which is distinct from the issue of 
symmetry restoration, is however desirable if possible [14]; and in §4 we also show how 
to incorporate it naturally within the LMA. The Friedel sum rule is also satisfied in the 
modified perturbation scheme [14], where the single self-energy (excluding the Hartree 

the 



(2) 

piece) is parametrized in the form E(co>) ~ v4E (cu)[l — BT, 
second order self-energy constructed from non-interacting Green functions containing a 
shifted chemical potential /to- The parameters A,B,fio are then found [14] by requiring 
that S(cl>) has the correct high-frequency behaviour (in 1/uj), that the atomic limit is 
recovered, and that the Friedel sum rule is satisfied. While simple and practicable, the 
approach is however too crude to handle e.g. the Kondo regime (as evident for example 

(2) 

in the fact that for the symmetric AIM where B = and A = 1, £(co>) — S (u) reduces 
to straight second order perturbation theory and hence produces a Kondo scale that is 
algebraically rather than exponentially small in U). 



,(2), 



- 1 : with S^ 2) M 



3.2. Mean field 

Before proceeding we reprise briefly elements of pure MF that will be required in the 
following sections (full details are given in [1]). The MF propagators Q a {uj) may be 
expressed as G a {u) = [uj + — ti a — A(u;)] _1 with e m = e; + S° = ei + ^(n — cr|yu|); and 
n, found self-consistently at pure MF level via (n^ ± n^) respectively (with (...) 
a MF average). The corresponding spectra D®(u) = D^{uo\e^x) are given (for the 
wide-band host explicitly) by 

„. . Anvr- 1 
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where x = |C7 and e; = e x + |n. The MF charge and moment are thus found from 
ion of 

^2 a / du D° a (u;e u x) (3.14a) 



2 

self-consistent solution of 



oc 



/■0 

2j / du D° a (u;e u x). (3.14b) 

„ J — oo 



We choose for later convenience to work at fixed U; in which case the MF \fi\ follows 
by solution of equation ( 3.14a|) for given ei, the charge n follows directly from equation 



( |3.14b|) , and the corresponding 'bare' e; then follows immediately from e; = e- x — ^n. 

The boundary to local moment formation at MF level is readily deduced from 
equations (|3.14a| ,b) and is given in closed from by 



2 

±- 

7T 



tan" 1 ( \IU' C 



U> 



UL 



(3.15) 



(which recovers figure 4 of [1]), where U' = U/tt and rj is the asymmetry (equation 
( p.2[) ); such that the MF > for any given < 1 when U' > U' c (or equivalently 
for all \rj\ < \r) c \ and any given U' > 1). 

4. LMA: practice 

In this section we consider first the approach to solving the problem generally within the 
LMA framework, for an essentially arbitrary approximation to the self-energies E^cu); 
before turning in §4.1 to the specific class of diagrams retained in practice, and the key 
consequences of stability and symmetry restoration. 

The LMA self-energies Ti a {uj) (equation (3.3)) are functionals of the MF 
propagators, viz Q a {uj) = [uj + — e ; + ax — A(co>)] -1 (with corresponding spectral density 
D®{u) = D^(w;ei,x), equation (3.13)). With interactions included beyond MF level, 
x = ^U\fi\ and e\ are determined self-consistently as below and naturally differ from 
their pure MF values; so for clarity the latter are denoted from now on as x = ^U\/iq\ 
and ef = e; + ^n (as obtained from the pure MF equations (3.14)). Equation (3.3) for 
T, a (uj) is given by 

Z a (u) = ^[n-*\ii\] + V ( ,(u) (4.1) 

where the first term is the static bubble diagram, E° = U f° du> D°_ a (uj; e;, x); and 
\fl\ = |/2(ei, x)\ and n = n(ei,x) are given by 

|/2| = ) a / du D° a (uj; e h x) (4.2a) 



|n| = ^/ du D" a (u;e h x). (4.2b) 
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It is of course the 'post-MF' contribution, T, a (u), that is all important; and a suitable, 
naturally approximate choice for which determines the extent to which the key physics 
of the problem is captured in practice. That is considered in §4, but for the present 
assume E cr (o;) = S cr [{^ cr (c<j)}] to be given. It depends generally on U (via the 
diagrammatic interaction vertices) and, since it is a functional of {Qa}, upon e\ and x; 
i.e. Tl a (uj) = T: a (u; ei, x) (with the U-dependence implicit). The symmetry restoration 
condition equation ( |3.8D is then of form: 



Ef(w = 0;e h x) - = 0;e is x) = U\jl(e h x)\. (4.3) 

And the Friedel sum rule equation ( |3.11] ) (equivalent to the Luttinger theorem [36]) is 
given by 

e i + jH^ux) - a\fl{ei, x)\) + = 0;e u x) = A tan[f (1 - n imp )] (4.4) 

(independently of a). 

The problem may then be solved in the following natural way, considered for (any) 
fixed U: 

(i) For given e; [or x], the symmetry restoration condition equation ( |4.3| ) is solved for 
x = \U\[i\ [or ej. 

(ii) n- imp is then obtained by solution of equations (3.12). 



(iii) All quantities in equation (|4.4j ) save e\ are then known, whence the bare e x follows 
directly from equation (|4.4|) . 

Alternatively, if one wishes to specify a bare e x from the outset, the procedure may 
be repeated by varying e Y until equation (Op is satisfied. The method above is 
straightforward, and in particular we remind the reader that (see §3) satisfaction 
of symmetry restoration (step (i)) guarantees the requisite low-energy Fermi liquid 
behaviour. As considered in practice (§4. Iff) it is also numerically rapid, with each 
of the iterative steps (i) and (ii) typically found to converge after a small number of 
iterations. 

Using the symmetry Q a {uj\e h x) = —Q_ a {—uj;—e^x) together with the general 
diagrammatic structure of T> a (uj) = it can moreover be shown that 

G a (u;ei,x) = — G- a (— u; — e\, x) and £ CT (co>; e- u x) = U — £_o-(— uj\ — e- u x)\ and in 
consequence that if (x, n imp , e{) is a solution of the above equations for ei = a, then 
(x, 2 — n imp , — [e; + U]) is also a solution for e± = —a. In other words the particle-hole 
transformation mentioned in §2, under which e\ — > — [e; + U] and n imp — > 2 — ni mp , 
corresponds to e- x — ej (with D a {uj;e\) = D_ a (— u; — [e; + U]) such that equation 
(2.11) follows directly, since D(u) = \ ^ a D a (uj)). From this it follows in turn that 
only e ; > need be considered (corresponding to n imp < 1); and that the particle-hole 
symmetric AIM, for which e\ = — y and n luip = 1, corresponds to e\ = 0. For the latter 
case, steps (ii) and (iii) above are redundant, the Friedel sum rule equation ( ]4.4| ) being 
guaranteed by particle-hole symmetry since S^(u; = 0) = | in that case; and solely the 
central symmetry restoration condition equation (|4.3j ) need be solved. The problem then 
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becomes that considered previously [27,28] (noting that E CT (o;) used therein is defined to 
exclude the Hartree contribution of y), with one minor exception; namely that in [27-32] 
\fx\ itself rather than \p,\ was employed on the right hand side of equation (|4.3j ). (The 
latter, used here, is technically correct; although we have confirmed that the results of 
[27-32] are entirely unaffected by this replacement.) 



4-1. Self- energies 



The specific class of diagrams contributing to £ CT (u;) = ^^[{^o-}] which we retain in 
practice is precisely that considered hitherto for the symmetric AIM [27-32]; shown in 
figure 1 and translating to 

du 1 



2vri 



(4.5) 



Physically, these diagrams capture dynamical spin-flip scattering processes: in which 
having, say, added a cr-spin electron to a — a-spin occupied impurity, the latter hops 
off the impurity and thus generates an on-site spin-flip (reflected in the transverse spin 
polarization propagator U~ cra (uj)) , before returning again at a later time. At the simplest 
level, likewise considered here, the polarization propagator is given by an RPA-like 
particle-hole ladder sum in the transverse spin channnel; specifically 

W- a (uj) = °n ff - ff (w)[l- U°W- a {u))- 1 (4.6) 

with °n cr ~ f7 (cG>) the bare particle- hole bubble, itself expressed in terms of the broken 
symmetry MF propagators {Q a }- Further details regarding the S CT (u;) diagrams and 
their physical content are given in [27,32] (see e.g. figure 3 of [32]). The full LMA 
self-energies Yi a (uj) follow directly from equations ( |4.1| ,5). 

Retention of the above diagrams for Y> a (uj) is moti- 
vated primarily on physical grounds, since they embody 
the dynamical coupling of single-particle excitations to 
low-energy transverse spin fluctuations that is essential to 
capture in particular the strong coupling Kondo regime. 
Other classes of diagrams are readily retained, notably 
those shown in figure 9 of [27] involving repeated particle- 
particle interactions in the transverse spin channel, as well 
as the repeated 'bubble' sum. But we find these to have 
little effect on the results given here, and that the class re- 
tained appears to handle rather well essentially all regimes 
of the AIM, from Kondo to empty orbital; reflecting at 
least in part the fact that the present approach asymptot- 
ically recovers straight second order perturbation theory 

in U in the weak coupling domain appropriate to the empty orbital regime (as will be 
shown explicitly in §5.2). Finally, although the atomic limit Ao = is not perturba- 
tively connected to the generic case of non-zero hybridization strength, we add that the 



Z a (co) = 




Figure 1. 

Principal contribution to the 
LMA £ CT (u;), see text. Wavy 
lines denote U. 
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present approach recovers it exactly. Here T, a {ui) = and simple MF alone is exact, 
encompassing both the n = (e; > = Ep) and n — 2 (ei + U < 0) regimes as well as 
the doubly degenerate local moment state arising in the singly occupied regime where 
e ; < < e ; + U . This of course is just a simple consequence of the fact that we work in 
general with broken symmetry MF propagators; but we note that recovery of the atomic 
limit is not at all a trivial matter if one seeks to obtain it via conventional perturbative 
methods, save for the accidental case of the symmetric AIM [2] where second order 
perturbation theory in U happens to be exact in the atomic limit. 

Since the diagrams for S„(w) have the same functional dependence on {Q a } as 
considered in [27] for the symmetric AIM, most of the specific analysis of §s 2.2, 3.2 of 
[27] goes through unaltered. We thus reprise only two required results, before focusing on 
the central issues of stability and symmetry restoration. First, independent of particle- 
hole symmetry, the (and hence 11°-° (u)) satisfy °n + ~(cj) = °II~ + (-a;). 
Only one polarization propagator need thus be considered, say °U + ~(uj); and for later 
use below, lm°U (T ~ CT (to) > V u, and vanishes linearly in \cu\ as cu — ► 0: 

-Im°U + -(u) |w|D?(0)£>2(0). (4.7) 

7T 

Second, the real/imaginary parts of °n + ~(ci;) and U + ~(uj) are related by the Hilbert 
transform 

= r MtM^M. (4 . 8) 

Using this, together with T\r aa (uji) = n <T_ °'(— u>i), equation ([4.5|) may be expressed as 
= U 2 Imn ff - CT ( Wl ) [eMGZA"! + + ei-untf+rfa + w)] (4.9) 



g^ )= l (4.10) 



where 



u — uji ± i0" 

denote the one-sided Hilbert transforms such that Q a {oS) = Q„ (uj) + Q~ (u); and 6{oj) is 
the unit step function. 

We turn now to the important issue of stability, namely (from equation fl4.8Q ) that 
7rRen +_ (a; = 0) = dui Imn +_ (co' 1 )/|ci; 1 | > of necessity. For this to be satisfied, 
using equation (|4.6|) and that Im°n +_ (c<j = 0) = (equation (|4.7|)), it follows that 

< URe°U + -(uo = 0) < 1 (4.11) 

is required. An explicit expression for Re°n+-(0) can however be obtained in direct 
parallel to that in [27] for the symmetric AIM; and is given by 

URe°U + -(u; = 0) = (4.12) 

with |/2(ei,x)| given by equation (4.2a) (and x = ^U\fi\). For (any) given U, the contour 
URe°U + -(uj = 0; e- x ,x) = 1 in the (e;, |/i|) plane encloses a region of instability where 
URe°U + - (u = 0; e^ x) > 1 and the positivity condition Ren +_ (c<j = 0) > is violated; 
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Figure 2. Dashed line: stability border in the (ej, plane for fixed U — 4ir. Solid 
line: corresponding solution of the symmetry restoration equation (4.3). Full discussion 
in text. 



outside this region by contrast, UKe U + ~ (u = 0; e^x) < 1 and stability is guaranteed. 
This is illustrated in figure 2 which, for U = U/A = An (and the wide-band AIM in 
practice), shows the stability border in the (e i? plane. 

Three important points should be noted here, (i) First, the divergences alluded to 
in §3 that bedevil standard diagrammatic resummations based on perturbation theory 
in U — and thus employing restricted HF propagators (where = is enforced) - 
arise from violation of the stability condition and hence analyticity; as is evident from 
figure 2 for e.g. the symmetric AIM (corresponding to e\ = as explained above). This 
deficiency, reflecting the associated instability of the restricted HF saddle point to local 
moment formation, does not arise in the present approach; which is why diagrammatic 
resummations of essentially standard form, but expressed in terms of broken symmetry 
MF propagators Q a , may be used with impunity, (ii) The stability border, illustrated in 
figure 2, corresponds precisely to the pure MF local moment, denoted by \fj,o\ = |/-fo( e i)|; 
for at pure MF level the local moment is obtained for any U as the self-consistent solution 
of \/jl \ = \fl(e i: |C/|// |)| (equation (3.14a)), and hence U°TI + ~(lu = 0;ei,x ) = 1 from 
equation (4.12). In consequence the transverse spin polarization propagator II + ~(c<j) 
(equation (4.6)) contains a pole at to — 0, arising physically because the pure MF state 
is a degenerate doublet with no energy cost for a local spin flip. This behaviour is 
correct for a local moment phase [32], and thus in the present context for the singly 
occupied regime in the atomic limit Ao = 0; but it is not of course correct for the 
ubiquitous Fermi liquid state characteristic of the metallic AIM with A ^ 0, where the 
characteristic energy for spin flips is on the order of the Kondo scale. The key point here 
however, is that this behaviour is specific solely to the pure MF level of self- consistency: 
from equation (4.12), URe°U + ~(uj = 0;ei,x) = 1 only if is determined by the pure 
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MF self-consistency equation, such that = \/j, (e{}\. (iii) And within the LMA it 
is the symmetry restoration condition equation (4.3) that, quite generally, determines 
the local moment (or equivalently x = ^U\/i\), as explained in the previous section. 
In consequence, ImIT + ~(c<j) will contain not a zero-frequency spin-flip pole but rather 
a resonance centred on a non-zero frequency cu m . This is the Kondo scale. Its origin 
within the LMA thus stems from self-consistent imposition of symmetry restoration; 
and since the latter guarantees Fermi liquid behaviour at low energies as explained 
in §3, its physical content is that it sets the timescale r ~ h/u m for restoration of the 
broken symmetry/degeneracy endemic at pure MF level. The solutions of the symmetry 
restoration equation (4.3) in the (e;, plane are illustrated in figure 2 for U = Air; and 
are also seen correctly to satisfy the stability criterion (as is always found in practice). 
We add further that the above behaviour of ImIT + ~(c<j), arising for the AIM generically, 
is illustrated e.g. in figure 2 of [27] for the symmetric model; and that the analytic 
description of Imll + ~(uj) given by equations (2.34,35) of [27] applies mutatis mutandis 
to the asymmetric case. 

Finally, referring to equation (4.9) for Yi a (u), it follows using \m.Qf(u)) = 
TnD° a (Lu)6(±uj) that 

p | uj I 

El(u) = e(-uj)U 2 / dc^x ImW-^uj^D^^ + u) 
Jo 

r° 

+ d(uo)U 2 Au 1 h^-' T {u 1 )D. <T {u 1 +u). (4.13) 

such that S^(cj) > ( since ImrP a {uj) and D®(u) are positive semidefinite), as required 
by analyticity. Using URe°U + ~ (uj = 0) < 1 (stability), it also follows from equations 
(4.6,7) that ImIT 7 " '(u) itself vanishes linearly in \uj\ as uj — > 0; and hence from equation 
(4.13) that S^cu) ~ 0{uj 2 ) as uj — > 0, as stated/used in §3 (equation (3.6b)). 

5. Results 

We now turn to results arising from the LMA in practice, i.e. with Yi a (uj) given 
approximately by equation (4.9). The strong coupling Kondo regime, including the 
question of universal spectral scaling therein, is considered in §5.1; and the mixed valence 
and empty orbital regimes in §5.2. The results are obtained straightforwardly: equations 
(4.3,4) and (3.12) are solved via the direct scheme described at the beginning of §4, the 
self-energies T, c {uj) then follow via equation (4.1), and the impurity Green function and 
hence single-particle spectrum D{oj) from equation (3.2). 

By way of overview figure 3 shows the resultant nA D(uj) vs uj = uj/ A on all energy 
scales, for = ei/A = —4 and U = U/A = 12, corresponding to an asymmetry 
(equation (2.2)) of rj — |. The strong coupling Kondo limit, where spectral scaling 
arises and the low-cu behaviour of the AIM maps onto the Kondo model, corresponds 
strictly to n; mp — > 1 (§2); although the Kondo regime is often regarded more loosely 
as spanning charges in the interval 0.8 < n imp < 1 (with n imp = n for the wide-band 
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Figure 3. Single-particle spectrum ttA D(uj) vs w = w/A for e; = e;/A = —4 and 
U = U/Aq = 12 (corresponding to an asymmetry r\ = |). The MF spectrum is shown 
for comparison (dotted line). Insets: low-energy Kondo resonance on progressively 
expanded scales. 



model considered). The example shown in figure 3 pertains to the Kondo regime, with 
m mp ~ 0.93; and the central Kondo resonance is shown on progressively expanded scales 
in the insets. The spectral asymmetry on essentially all energy scales is evident: in the 
expected relative disposition of the high-energy Hubbard satellites centred on oo ~ — | e il 
and U — |e;|; and in the overall shape of the Kondo resonance (left inset), including the 
lowest scales (right inset). In the latter case, since the Friedel sum rule is fully satisfied, 
tiAqD(uj = 0) = sin 2 [|ni mp ] is correctly recovered; and since rii mp < 1 the spectral 
maximum in D(u>) is then pushed slightly above the Fermi level, in accordance with the 
quasiparticle form equation (2.7). 

The qualitative deficiencies of the pure MF spectrum (also shown in figure 2) are 
directly apparent, even on the high-energy scales of the Hubbard satellites. These show 
clearly the effects of additional many-body broadening over and above the purely elastic 
scattering processes captured at MF level; and the simple physical origin of which is the 
same as that described in [27] for the symmetric AIM, leading as seen in figure 3 to a 
width doubling (and peak intensity halving) of the satellites relative to pure MF. 

5.1. Kondo scaling regime 

Figure 4 illustrates the strong coupling Kondo regime where 7ii mp — > 1. In figure 4a 
spectral evolution is shown at a fixed asymmetry 77 = l—2\t]\/U = 0.2; with progressively 
increasing U = 25, 30, 35 such that both Hubbard satellites move steadily outwards 
with a fixed ratio of their peak maxima. In figure 4b by contrast, 1\ = —10 is fixed and 
U = 2|e;|/(l — 77) is increased from the symmetric limit rj = through 77 = 0.1, 0.2, 0.3 
and 0.4; such that the lower satellite is 'frozen' while the upper, centred onw~(/- |ei|, 
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Figure 4. (a) ttA D(lu) vsii = uj/A for fixed asymmetry r\ = 1 — 2|?i|/f7 = 0.2 and 
progressively increasing U — 25, 30, 35. (b) Spectral evolution for fixed e; = —10 and 
progressively increasing U — 2|e;|/(l — rj) from rj = (symmetric limit) through 77 = 
0.1, 0.2, 0.3, 0.4. Insets in either case: Kondo resonances vs Co, with increasing {7's 
from outside to inside. 

moves to progressively higher energies. In both cases shown however, since ni mp = 1 for 
all practical purposes, D(u) is peaked at the Fermi level uo = such that nA D(0) = 1, 
as evident from the insets which show the central Kondo resonance on the half-height 
scale. And with progressively increasing U in either case the Kondo scale, reflected in 
the width of the resonance, is rapidly decreasing. 

The first question then is the behaviour of the Kondo scale in strong coupling, 
equivalently the spin-flip scale u m determined (see §4) by symmetry restoration. This 
may be obtained analytically within the present LMA; using initially, as now outlined, 
a generalization of the arguments detailed in [27] for the symmetric model, and with 
reference to equation (4.9) for E <7 (a;). In the strong coupling Kondo regime, the spectral 
weight of ImII + ~(cc>) is confined entirely to w > 0, and J °°(dc<j/7r) ImII + ~(c<j) = 1; 
which behaviour reflects physically the saturation of the local moment — > 1). The 
resonance in Imll + ~(w) is centred on the low-energy Kondo spin-flip scale u> m , and on 
scales of this order G^(ijS) is slowly varying; whence equation (4.9) for Y^(uo = 0) reduces 
asymptotically to S^(0) ~ U 2 ReQ^(u m ). But ReQ^(uj) is given by the one-sided Hilbert 
transform equation (4.10), which as uj — > is dominated by the log singularity arising 
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necessarily because the host is metallic {D g Auj = 0) 7^ 0). This uj — > asymptotic 
behaviour is captured by ReQ7(u>) ~ Z)J(O) ln[A/|a;|]; where a high-energy cutoff A of 
order min[[7, D] is employed (D here being the host bandwidth, if one wishes to retain 
it as finite). The precise value of the cutoff is immaterial in the following analysis, the 
key point being that the prefactor to the log divergence is precisely D®(0). And in 
strong coupling (where x = hU\fj,\ 



(3.13) by 7rD{(0) 
£f (u = 0) is: 



-> \U), -D°(0) is given asymptotically from equation 
Hence the strong coupling asymptotic behaviour of 



0) 



u 2 



In 



vr (e, + \UY 

A directly analogous argument may be used for Xf(w = 0), leading to £f(0) 
U 2 ReQ^(—uj m ) ~ — U 2 D®(ti) In [A/u; m ]; and hence via equation (3.13) to: 



(5.1a) 



£> 



0) 



f/ 2 



In 



A 



(5.1b) 



7T ( ei - \U) 2 |W n 

The C/- and ej-dependence of the Kondo scale spin-flip scale u m now follows from 
the symmetry restoration condition equation (4.3), which in strong coupling (where 
\p,\ — > 1) reduces simply to S^(0) — Sf(0) = U; and hence using equation (5.1): 

. ( ~* [ef-jf\ .... 

^~ Xe ^{2^u^Tf\)- (5 ' 2) 

For the symmetric AIM (77 = 1 — 2\ei\/U = 0), where e; = (§4), equation (5.2) alone 
yields the Kondo scale; viz u m ~ A exp[— nU/8A ] where the exponent is exact [27]. 
To determine u m in general, the dependence of e; on the bare parameters e; and U is 
naturally required. This may be obtained using the Friedel sum rule equation (4.4) 
which, in the strong coupling regime where n imp , n and \p,\ — > 1, gives the asymptotic 
behaviour ^f(u> = 0) ~ |ei| — U. Combining this with equations (5. lb, 2) yields a simple 
quadratic for e; with solution e\ = [U — 2|ei|]/(4/); where / = f(rf ) is given by 



'-5 



1 + (1 



v 2 ) 1 



(5.3) 



such that / G [1, |] for the i] G [0, 1] domain relevant (§2) to the Kondo regime. The 
resultant is now used in equation (5.2) for the Kondo scale u m , to give 

7r|ej|(£/-H) 1 
2A U f, 



0O n 



A exp 



(5.4a) 



A exp 



(5.4b) 



2pJ/) 

where equation (2.13a) for pj appropriate to the Kondo model is used. 

Equation (5.4) gives the strong coupling Kondo scale for the AIM arising from 
the present LMA (as we have confirmed numerically). The exponent therein differs in 
general, by the factor of f(i] 2 ) G [1, |], from the exact result for the Kondo model; 
which is given to leading order in pj by [2] u m oc exp(— l/2pj) independently of the 
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Figure 5. Kondo scaling regime, (a) ttA D(uj) vs u)/u m for fixed asymmetry 
r] = 1 — 2|ei|/f/ = 0.2 and increasing {/ = 25, 30 , 35 (cf figure 4(a)). The spectra 
collapse to a common scaling form, (b) Resultant scaling spectrum for r\ = 0.5 (solid 
line); also shown are the \ui\/oj m > 1 tails of the scaling spectra for r) = (dotted line) 
and r\ = 0.2 (dashed line). Inset: scaled Kondo resonances at higher resolution for 
7] = (dotted line), 0.2 (dashed line) and 0.5 (solid line). 

strength of potential scattering embodied in Kj J = rj. As such it is exact only for the 
symmetric model where 77 = and / = 1 (although we add that / is slowly varying in 
77, lying e.g. within 10% of unity for 77 < 0.6). Nonetheless we regard recovery of an 
exponentially small Kondo scale, approximate in general but close to the exact result 
in an obvious sense, as a non-trivial outcome of the present theory. Moreover, provided 
the resultant Kondo scale is exponentially small in strong coupling (so that a clear 
separation between high and low energy scales arises), its precise dependence on the 
bare material parameters is nigh on irrelevant to the central question of scaling in the 
Kondo regime; viz the behaviour of the single-particle spectrum D{u) = F(uo/uo m ) as a 
function of oo/uj m in the formal Kondo limit u ni — > 0. This is known from previous work 
on the symmetric AIM [28]. Here, although the exponent in equation (5.4) for u m is 
exact, the prefactor A is simply a high-energy cutoff and certainly approximate; and yet 
in comparison to NRG calculations [33] the resultant universal scaling spectrum of the 
symmetric model is captured quantitatively by the LMA (see e.g. figures 2,3 of [28]). It 
is to the questions of scaling for the asymmetric AIM that we now turn. 
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Spectral scaling 

The first issue in regard to spectral scaling in the Kondo regime is: does the 
asymmetry/potential scattering rj = K/J play a role? At first sight one might be 
tempted to answer no; arguing, pace for example the Kondo Hamiltonian equation 
(2.12), that potential scattering can be largely eliminated via a suitable canonical 
transformation of the host band {cj^}. And as far as the scaling spectrum is concerned 
this is correct at the Fermi level uj = 0; where, as explained in §2 and in accordance with 
the quasiparticle form equation (2.7) for G(uj), ttAqD(uj = 0) = 1 provided n; mp = 1 but 
regardless of how the Kondo limit rii mp — > 1 is reached (i.e. independently of rj). Indeed if 
one extrapolates the quasiparticle form equation (2.7) beyond its domain of applicability 
as the limiting low-a; behaviour of G(u), as e.g. in microscopic Fermi liquid theory to 
lowest order [2], the resultant scaling spectrum D(uj) = F(uj/ A Z) for n imp = 1 (where 
e[ = 0) is a simple Lorentzian centred on uj = and with no 77-dependence; which 
conclusion arises also if D(uj) is approximated by the local spinon spectrum [25]. 

If however the effects of potential scattering were entirely irrelevant, then the scaling 
spectrum would obviously coincide with that for the symmetric model (77 = 0) [33], and 
as such D(uj) = D(—uj) would be symmetric about the Fermi level for all uj/uj m (or 
uj/ A Z) and any r\ G [0, 1]. Such behaviour is in our view rather implausible, and we 
are not aware of a convincing general argument for it. Within the usual U = 00 NCA [17- 
19] with n imp ~ 1 for example, the low-energy Kondo resonance is certainly asymmetric 
in uj (albeit not correctly centred on uj = as n imp — > 1 since the Friedel sum rule is 
not consistently satisfied [38]). Neither would it square readily with NRG results [39] 
for the asymmetric AIM, which show an evident asymmetry in the tu-dependence of the 
low-energy D(uj) in the Kondo regime, see e.g. figure 5 of [39]; the only counterargument 
to which being (we believe improbably) that the latter calculations, obtained for fixed 
U = An by varying ej = — |ej|, are insufficiently close to the Kondo limit to eliminate 
asymmetry in the Kondo resonance. 

If then the potential scattering/asymmetry embodied in rj = 1 — 2\ii\/U = K/J is 
not irrelevant, the obvious question is how it influences the scaling behaviour of D(uj). 
Within the LMA as now shown, the answer is that a continuous family of universal 
scaling spectra arise; D(u) = F(uj/uj m ) exhibiting one-parameter scaling in terms of 
oo/uj ni for each rj e [0, 1] (which is consistent with the known line of RG fixed points [9], 
one for each K). For fixed 77 = 1 — 2|ej|/£/ = 0.2, figure 4a above shows ttAqD(uj) vs the 
'absolute' frequency uj = uj/ A , for increasing U = 25, 30, 35; with the corresponding 
behaviour of the Kondo resonances shown in the inset. As shown in figure 5a, the Kondo 
resonances for the same data collapse to a common form when expressed in terms of 
uj/uj m , yielding the i] = 0.2 scaling spectrum discussed further below [we add here only 
that, as discussed in [27,28], the small feature at uj/uj m ~ 1 is entirely an artefact of 
the specific RPA-like form for the ImII +_ (w) employed here; it can be removed, but we 
are content to live with it in the following]. Collapse to a scaling form occurs for each 
rj with increasing U, and figure 5b shows the 77 = 0.5 scaling spectrum (solid line). The 
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inset thereto compares the central portion of the scaled Kondo resonance for rj — 0, 0.2 
and 0.5; from which the increasing asymmetry with rj is evident, modest though it is 
on the scales shown. The asymmetry induced by increasing rj is more clearly evident in 
the 'tails' of the scaling spectra, which for |a>|/u; m > 1 are also shown in figure 5b for rj 
= 0, 0.2 and 0.5. 

The first point to note about the LMA scaling spectra is that they are not of the 
Lorentzian form suggested by simplistic approaches. While the quasiparticle form is 
correctly recovered at sufficiently low frequencies |u;|/u; m <C 1 (§3), the scaling spectra 
are by contrast dominated for |u;|/u; m > 1 by the long, slowly varying tails evident in 
figure 5. From a recent LMA study of the symmetric AIM in strong coupling [28], these 
are known to exhibit a very slow logarithmic decay, and to give excellent agreement 
with NRG scaling spectra [33] (see e.g. figure 2 of [28]). 

The logarithmic tails naturally persist when rj ^ 0, and in parallel to [28] it 
is likewise possible to determine analytically their \uj'\ = \cu\/u> m ^> 1 asymptotic 
behaviour. For uj' = \uj'\ 3> 1 we thereby find 



the (asymmetry) ^-dependence of which is embodied in the factors g±(rj) = f (i] 2 )[l±ri]~ 1 
(with / from equation (5.3)). For negative u and \uj'\ ^> 1 by contrast, 7rA -D(u;) is 
again given by the form equation (5.5), but with g + and g_ simply interchanged. In 
the symmetric case where g±(0) = 1, equation (5.5) reduces to the result compared to 
NRG calculations in [28]. It reproduces quantitatively the uj'- and 77-dependences of the 
spectral tails shown in figure 5 for \uj'\ > 10 or so, and is quite satisfactory down to 
\lu'\ ~ 2; although since 77 = 1 — 2|e i |/£7 = K/J involves a ratio of 'bare' parameters, 
the ^-dependence of g± is undoubtedly approximate save for the symmetric case. 

Finally, we add that since universal spectral scaling arises for fixed 77 = 1 — 2\i[\/U 
upon increasing U, it does not therefore arise on increasing U for fixed |ei| = — e ; (or vice 
versa). For example the Kondo resonances illustrated in figure 4b for fixed |ei| = 10 and 
increasing U do not collapse to a common scaling form when expressed as a function 
of uj/u m , since each corresponds to a different asymmetry rj. We naturally expect 
the same conclusion to apply at finite T, namely that D(uo\T) = F(uj/uj m ;T/uj m ) will 
exhibit universal scaling only for fixed rj; and hence also to the T/u m scaling behaviour 
of transport properties such as the resistivity, since these are determined (see e.g. [39]) 
via transport integrals in which the transport rate r _1 (c<j; T) oc D(u; T). 

5.2. Mixed valence and empty orbital regimes 

Although we have focussed latterly on scaling behaviour in the Kondo limit, the Kondo 
regime of the AIM is usually regarded more generally as encompassing charges in the 
interval 0.8 < ni mp < 1; while the mixed valence (MV) and empty orbital (EO) regimes 
of behaviour correspond respectively to 0.3 < n hnp < 0.8 and n imp < 0.3. Here we 
consider spectral evolution on all energy scales upon traversal of the Kondo, MV and 



itA q D(uj) ~ - 




g«7 + ln|u/|] 2 + l [^_lnM] 2 + [l + 4<?_P 



7T 



1 (l + 4<7_) 
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Figure 6. Kondo regime. ttAoD(uj) vs lj = iv/Ao for fixed U — 4tt and e; = — ^ 
(dotted line), -4, -3 and -2. Inset: corresponding Kondo resonances vs cD. 

EO regimes by progressively increasing e; from the symmetric limit e; = —U/2; for a 
fixed interaction strength of U = An, as considered in the NRG calculations of [39]. 

The Kondo regime is illustrated in figure 6 where 7iA D(u) vs uj = u/A is shown 
for e ; = — U /2, -4, -3 and -2 (with n imp ~ 0.8 in the latter case and the resultant charges 
n imp found to agree well with NRG results [39,38], to which they are compared in figure 
7 below). We regard the agreement between figure 6 from the present LMA, and the 
corresponding NRG results in figure 5 of [39], as rather good. The principal differences 
between the two reside first in the Hubbard satellites which, we believe correctly, are 
more pronounced in the LMA; it being known that NRG does not do full justice to the 
many-body broadening of the Hubbard satellites, rendering them too diffuse [40]. And 
second, consistent with the expectation above, on the low-energy scale of the Kondo 
resonance (figure 6 inset) the LMA does not capture quantitatively the absolute value 
of the Kondo scale reflected e.g. in the resonance width; although the relative evolution 
of the Kondo resonances with increasing ii is respectably captured. 

On further increasing e; the MV regime is entered, where the impurity charge 
n = n imp drops quite rapidly with increasing e ; as the EO regime (ei > 1-2) is 
approached; as evident from figure 7 where the LMA n imp vs e; is shown, and compared to 
corresponding NRG results [38,39]. Spectral evolution in the MV and into the EO regime 
is illustrated in figure 8 for e ; = —1, 0, +1, +2, and likewise exhibits the characteristic 
behaviour found in NRG calculations (see e.g. figure 8 of [39]). The Kondo regime is of 
course typified by a clear separation of energy scales, reflected in the exponentially small 
quasiparticle weight Z oc u m / A Q and hence exponentially narrow Kondo resonance. This 
behaviour is lost quite rapidly on entering the MV regime, where Z rises to become of 
order unity and the width of the Kondo resonance correspondingly becomes of order 
A , as evident from figure 8. This is accompanied in turn by 'loss' of the lower Hubbard 
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Figure 7. Local impurity charge n = n; mp vs e; for fixed U = A-k. Solid circles: LMA. 
The solid line is a guide to the eye. Open squares: NRG results (from [39] for e; < 1 
and [38] for e; > f). The non-interacting limit charge is also shown (dotted line). 

satellite, which becomes a barely perceptible low-energy shoulder; and by concomitant 
intensity erosion of the upper Hubbard satellite, although the latter remains a high 
energy feature centred on ~ e; + U . 

While the symmetry restoration condition (equation (4.3)) fundamental to the LMA 
is always found to be satisfied, we note that in the lower-n portion of the MV regime and 
into the beginning of the EO regime we encounter difficulties in solving the Friedel sum 
rule fully self-consistently (as in the algorithm specified in §4). This occurs in practice 
when |E^(0)| is small (< 0.3); and while the reasons for it are largely technical we do 
not find it surprising since in this non-universal regime comparably small contributions 
to S^-(O) (in addition to the spin- fluctuation diagrams retained) will undoubtedly play a 
role. In practice, as relevant to the = 0-2 examples in figure 8, we have circumvented 
the matter, simply by replacing n imp in equation (4.4) (or equation (3.12a)) by its MF 
counterpart n = n(e;,x); and although the Friedel sum rule is not then fully satisfied, 
the differences between n and n irap determined by spectral integration (via equation 
(3.12b)) are modest (e.g. ~ 10% for ej = +1). 

The EO regime itself also illustrates an important facet of the LMA, for although 
the approach is naturally designed to capture the strong coupling physics inherent to 
the Kondo regime, it is nonetheless perturbatively correct in weak coupling as the non- 
interacting limit is approached. This arises for fixed U on progressively emptying the 
impurity site by increasing ei further in the EO regime (where we add that the Friedel 
sum rule is in general satisfied to full self-consistency); as evident in figure 7 where 
nimp = n asymptotically approaches the non-interacting charge, shown as a dotted line. 
Corresponding spectral evolution, ttAqD(uj) vs uj = w/A , is shown in figure 9 for e; = 4, 
7 and 10. For e; = 4 the upper satellite remains apparent as a weak shoulder that all 
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Figure 8. Mixed valence regime. ttAoD(lj) vs Cj = oj/Aq for fixed U = Atk and 
h = -f, 0, +1, +2. 
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Figure 9. Empty orbital regime. itA D(lj) vs Cj — lj/A for fixed U = 4ir and 
?i = +4, 7, 10 (solid lines). For q = 10 the non-interacting limit spectrum is shown 
(dotted line). Comparison is also made in all cases to second-order perturbation 
theory (SOPT, dashed line); for e; = 7 and 10 the LMA and SOPT spectra are 
indistinguishable. 



but vanishes by e; = 7, and is entirely gone by i{ = 10 where the resultant spectrum is 
seen to be near coincident with that for the non-interacting limit (dotted line). 

In particular, the LMA recovers with increasing e\ the weak coupling result of 
second-order perturbation theory (SOPT) in U . The latter corresponds simply to 
enforcing x = ^U\fx\ = (so that So-(c<j) = S(a;), Q a (uo) = Q{u) etc. ), replacing Y\r cnT in 
equation (4.5) by the bare polarization bubble °n, and solving via steps (ii) and (iii) of 
the algorithm specified at the beginning of §4 (such that the Friedel sum rule is satisfied 
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to full self-consistency). For the LMA itself, x = jU\/j,\ is of course determined (step 
(i)) via symmetry restoration equation (4.3). The resultant local moment \fj,\ is thereby 
found to diminish with increasing e; in the EO regime and vanishes for e; = ej ~ 4.3 
for the chosen U — Att (far above the ej ~ at which, from equation (3.15), moments 
vanish at pure MF level; thus ensuring, as found for the symmetric AIM [27], that 
physical properties evolve smoothly as i{ passes through e io ). For e; > e io where = 0, 
the transverse spin polarization propagator equation (4.6) entering the LMA self-energy 
equation (4.5) may then be expanded perturbatively in U; the leading term of which is 
precisely the SOPT result for the self-energy. In figure 9 LMA and SOPT single particle 
spectra are compared for e; = 4, 7, 10; the SOPT spectra being given by the dashed 
lines. For e; = 4 the two are evidently very close (despite the persistence of moments 
within the LMA). And for e; = 7, 10 the LMA and SOPT spectra are to all intents and 
purposes coincident: both are shown in the figure, but the differences are imperceptible. 

6. Summary 

We have considered in this paper single-particle dynamics of the Anderson impurity 
model, a longstanding issue [1,2] pursued here via development of the local moment 
approach [27-32] to encompass the general asymmetric case. The LMA can handle 
interaction strengths from strong to weak coupling, as well as dynamics on all energy 
scales. And as implemented in practice, it appears to capture rather well the wide 
spectrum of physical behaviour inherent to the problem; from the non-perturbative 
Kondo physics characteristic of strong coupling, through mixed valence to the ultimately 
perturbative empty orbital regime. This in turn has ramifications beyond impurity 
models per se, to lattice-based systems within DMFT [5-8]; for the simple, practicable 
approach developed here can be extended further to encompass important problems such 
as e.g. the periodic Anderson model away from half-filling, and hence heavy fermion 
physics. 
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